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ABSTRACT 



Aims. We follow the evolution of galaxy systems in numerical simulations. Our goal is to understand the role of density perturbations 
on various scales in the formation and evolution of the cosmic web. 

Methods. We perform numerical simulations with the full power spectrum of perturbations, and with a spectrum cut at long wave- 
lengths. In addition, we have one model, where we cut the intermediate waves. We analyse the density field and study void sizes and 
density field clusters in diff'erent models. 

Results. Our analysis shows that the fine structure (groups and clusters of galaxies) are created by small-scale density perturbations 
of scale < 8 Mpc. Filaments of galaxies and clusters are created by perturbations of intermediate scale from ~ 8 to ~ 32 Mpc, 
and superclusters of galaxies by larger perturbations. 

Conclusions. We conclude that the texture of the cosmic web is determined by density perturbations of the scales up to ~ \00h~^ Mpc. 
Larger perturbations do not change the texture of the web, but modulate the richness of galaxy systems, and make voids emptier. 



Key words, large-scale structure of Universe - 
matter - Methods: numerical 
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1. Introduction 

The first studies of the three-dimensional distribution of 
galaxies and clusters of galaxies demonstrated that they are 
not distributed randomly, but are mostly located in filamen- 
tary superclusters connected by galaxy filaments to form a 



ies (lGre2orv & ThomDsonl 119781: iJoeveer & Einastol 


1978; 


'Joeveer et al.',' 


1978:"Tarenghi et al.Vl978VTullv & Fisheil 


1978; 


Einasto et al., 


1980; ZeldovichetaL, 1982). This observational 



picture was rath er similar to the theoretical prediction of 

■ IZeldovichl(ll970h . the s o-called pancake s cenario of galaxy for- 

■ mation, as discussed bv lZeldovichI (119781) . However, there were 
" some important difiTerences between the models and the ob- 
servations. In the real world, the most common structural ele- 

■ ments were filani ents of galaxies and their clusters, whereas the 
! lZeldovichl(ll978h scenario predicted the format ion of flat form- 
less pa nca ke-like systems (for a discussion see IZeldovich et al.l 
(1982) and Einas to et al.l (11983 )). The simulations performed for 
the neutrino-d ominated universe, called the Hot Dark Matter 
(HDM) model (Iwhite et al.Lll983l) . had similar problems. 

The problem of the absence of the fine structure was solved 
w hen the Cold Dark Matter (CDM) scenario was suggested 
by iBond et al.l ([1982') and others. The basic diflTerence between 
the HDM and CDM scenarios is the absence or presence of 
small-scale density fluctuations. The comparison of simulations 
for both dark matter types demonstrated the advantages of the 
CDM model - in this model, the fine structure is present in the 
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form of filanients o f various scales (see lMelott et al.l (Il983l) and 
IWhite et al.1 (Il987h among others). An even closer agreement 
of simulations with observations is obtained in mod els with a 
cosmological constant, ACDM, as demonstrated by iGramannI 
(1988). 

One problem remained: both numerical simulations and di- 
rect observations suggested that the dominant structural ele- 
ments of the cosmic web are fi laments of various scales, whereas 
the classical IZeldovichl(ll970h scenario predicted the formation 
of flat pan cake-like systems. This problem was solved when 
iBond et al] (|T996) studied in detail the formation of the large- 
scale structure. They found that primordial tidal fields play an 
important role in the evolution of structure, thus the dominant 
elements are filaments, not flat pancakes. 

iBond et aD (Il996h showed why the cosmic web has a fila- 
mentary character, but could not explain, why there is a large 
variety of systems from small groups and weak filaments to rich 
clusters and superclusters. It is clear intuitively that the main 
reason for the formation of systems of various richnesses is the 
presence of density perturbations of diflTerent scales. This idea is 
justified by comparing the HDM and CDM scenarios - the fine 
structure of the cosmic web is generated by density fluctuations 
of medium and small scales that are absent in the HDM model. 

The role of very large density waves on the for mation of the 
structu r e has been studi e d by a num ber of authors. iKnebe et al.l 
('2000'),"BaglaAPrasad| (|2006l), and'Bagla et all (|2009|) investi- 
gated the eflfect of the resolution of cos mological simulations 
on the properties of the filamentary web. Power & Knebd (l2QQ6h 
studied the impact of the simulation box size on the properties 
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Fig. 1. The luminosity density field of the SDSS in a spherical shell of 10 Mpc thickness at a distance of 240 Mpc. To 
enhance the faint filaments in voids between the superclusters, the density scale is logarithmic, in units of the mean luminosity 
density for the whole DR7. The rich complex in the lov ver area of the pictur e is part of the Sloan Great Wall; it consists of two very 
rich superclusters, SCL 111 and SCL 126 in the list by lEinasto et alJ (1200 ll) . The coordinates x, y are defined by eqs. (1) and (2). A 
cellular distribution of superclusters and filaments of various richness can be clearly seen. 



of the simulated cosmic web. ^ Bernardeau et al.l (l2002h gave a 
review of how non-linear perturbation theory helps us to under- 
stand the large-scale structure of the Universe and various statis- 
tical tools used to characterise its properties. 

It is well known that quasi-relaxed systems of galaxies, such 
as groups and clusters, are formed by perturbations, which have 
entered the non-linear stage of their evolution. The scale of the 
non-linearity is usually taken equal to about 8 Mpc. The rich- 
est non-percolating systems of galaxies are superclusters. Rich 
superclusters form a cellular distribution, with large voids sur- 
rounded by rich superclusters. The characteristi c diameter of 
these supervo ids is of the order of 100 Mpc (lEinasto et al.L 
Il994ll997bh . 

Supervoids are not empty but contai n a hierarchy of voids , 
as demo nstrated by many autho rs (iMartel & ^yassermanL 
fT993; "Lindner et al.', 1995; "Gottlober et al.", ^2003"; 
lAragon-Calvo et al., 2007; von Benda-Beckmann & Miiller, 



2QQ8I: Ivan de W evgaert et all \200^ Ivan de Wevgaert & Platenl. 
20091: lAragon-Calvo et al.L 



2010al: Ijones et a l.l '2010'). It 
is natural to expect that the web of rich superclusters is 
formed by density perturbations of scale 100 Mpc and 
above. The skeleton of the cosmic web is discussed by, 
among others, [ Sheth & van de Wevgaert (2004), Hahn et al. 
(120071) . iForero-Ro mero et al.. (.2009.) . Sousbie et al. ( 2008, 
120091). lArag on-Cal vo et al.1 (l2010bh . Bond et al. ( 2010b a), 
and lEinasto et al.i iQOlO). Thus, we can say that the role of 
small-scale and large-scale density perturbations is relatively 
well understood. However, it remains unclear how perturbations 
of intermediate scale combine to form elements of the cosmic 
web. 

We aim to investigate in greater detail than before the role 
of density perturbations of diff'erent scales in the formation of 



the cosmic web. We perform numerical simulations of the for- 
mation and evolution of the large-scale structure, using diff'erent 
box sizes and resolutions. Furthermore, to elucidate the influence 
of perturbations of difl'erent scales, a varying large-scale cutofl* 
is introduced into the power spectrum of initial perturbations. 
All simulations of a given series have identical initial conditions 
with random initial positions and velocities of test particles, but 
the amplitude of all perturbations on a scale exceeding a given 
one is forced to be zero. In this way, we can follow how systems 
of galaxies grow under the influence of perturbations of various 
scales. 

To characterise the efl'ect of varying the cutofl' scale of den- 
sity perturbations, we consider three quantitative tests: the mean 
radii of voids defined by groups and clusters of galaxies, the 
mass functions of clusters of galaxies, and the density distribu- 
tions of particles in the void and the supercluster core regions. 

The paper is organised as follows. In the next section, we 
describe numerical simulations used to model the evolution of 
the cosmic web. For comparison, we also calculate the SDSS lu- 
minosity density field. In Section 3, we analyse the cluster mass 
functions and void sizes. We analyse the density distributions in 
the void and the supercluster core regions and discuss our re- 
sults in Section 4. The last section presents the summary of the 
analysis. 



2. Modelling the evolution of the cosmic web 

2.1. The SDSS luminosity density field 

To compare simulated void and cluster data with actual data, we 
used the recently completed Sloan Digital Sky Survey (SDSS) 
Data Release 7 (DR7) (see lAbazaiian et"ai] (l2009i) ). In the 
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present analysis, we used only the contiguous northern zone 
of the DR7. The analysis was made in several steps. First we 
calculated the luminosity density field. We estimated the total 
luminosities of groups and isolated galaxies in a flux-limited 
sample, using luminosity weights that take into account galax- 
ies and galaxy groups too faint to fall into the observational 
window of absolute magnitudes at the distance of the galaxy . 
For details of the data r eduction, we refer to Ta go et all (120101) 
and lTempel et al.l (l201ll) . The high-resolution luminosity density 
field was calculated with the B3 spline of kernel size 1 Mpc. 

The high-resolution luminosity density field was found us- 
ing spherical coordinates (the SDSS coordinates rj and A) and 
distance. To obtain an impression of the luminosity density field 
of the whole northern region we show in Fig. [T] the luminosity 
density field in a spherical shell at a distance of 240 Mpc 
from ourselves. It is based on the SDSS angular coordinates t] 
and A, and the comoving distances Jgai of galaxies 

X = -dg^\A, (1) 
y = dgs,\r] cos A, (2) 
z = Jgai - do, (3) 

where do = 50 Mpc is the minimal distance used in the cal- 
culation of the density field. A supercluster catalogue based on 
the luminosity density field of the fu ll contiguous northern SDSS 
region was published by Liivamag i et al.l (|201Q|) . The represen- 
tation in spherical coordinates serves basically for illustrative 
purpose, as it is practically free of distance dependent selec- 
tion eff'ects. The spherical shell has a thickness of 10 Mpc 
to enhance the filamentary distribution of galaxies, clusters, and 
superclusters. The mean distance of the shell corresponds to 
the distance of the Sloan Great Wall, seen in the lower area 
of the figure. We see that it consists of two very rich super- 
clusters: SCLlll a nd SCL126 according to the catalogue by 
lEinasto et all (l200lh . 

2.2. Simulation of the cosmic web for density perturbations 
of various scales 

To understand the formation of the filamentary supercluster- void 
web correctly, we must perform numerical simulations for a box 
that contains both small and large waves. The smallest units of 
this network are galaxies, the most frequent systems of galax- 
ies are groups and clusters of galaxies. The characteristic scale 
of groups is I Mpc (galaxies are still about 10 to 100 times 
smaller), thus the simulation must have a resolution of at least 
this scale. On the other hand, the largest non-percolating systems 
of galaxies are s uperclusters, which have a characteristic scale 
of 100 h-^ Mpc (Port "1983'; 'Zucca et al., 1993; Einasto et al.", 
19941: iKalinkoy & Kuneva. 1995; Einasto et al, 1997a, 2001; 
Erdogdu et al.1. |2QQ4 lEinasto et al.L |2007|). Superclusters have 
rather diff'erent richness, from small systems like the Local 
Supercluster to very rich systems such as the Shapley 
Supercluster. It is clear that this variety has its origin in den- 
sity perturbations of still larger scales. Thus, to understand the 
supercluster- void phenomenon correctly, the influence of very 
large density perturbations should also be studied. 

To derive both a high spatial resolution and the presence 
of density perturbations on a large scale interval, we used a 
number of simulations in boxes of sizes from 100 Mpc 
to 768 Mpc, and various resolutions of A^^^.^ = 256^ and 

^grid ~ ^^^^ particles and simulation cells. To see the efl'ect 
of perturbations of various scales, we use simulations with the 



full power spectrum, as well as with a power spectrum truncated 
at wave-numbers kt, so that the amplitude of the power spec- 
trum on large scales is zero i.e., P(k) = 0, if ^ < kt, wavelength 
At = 271 /kt. 



Table 1. Parameters of the models. 



Model 


L 


^cut 


A^ciA 


A^ciB 




(1) 


(2) 


(3) 


(4) 


M768.768 


768 


768 


102783 




M768.128 


768 


128 


105496 




M768.032 


768 


32 


184193 




M768.012 


768 


12 


213569 




M256.256 


256 


256 


54224 


5578 


M256.064 


256 


64 


56000 


6077 


M256.032 


256 


32 


60985 


6996 


M256.016 


256 


16 


70813 


7624 


M256.008 


256 


8 




125 


M256.864 


256 


8-64 






L256.256 


256 


256 


341758 


59100 


L256.128 


256 


128 


343706 


59313 


L256.064 


256 


64 


351198 


61831 


L256.032 


256 


32 


372298 


68221 


L256.016 


256 


16 


426699 


84691 


L256.008 


256 


8 


546414 


126832 


LlOO.lOO 


100 


100 


54348 


173655 


L100.032 


100 


32 


57878 


182675 


L100.016 


100 


16 


62396 


197053 


LI 00.008 


100 


8 


71726 


237499 



Notes. 

Column 1: L - the size of the simulation box in hr^ Mpc; 
Column 2: /Icut - the cut-off scale in h'^ Mpc; 

Column 3: The number of DF clusters in the high-resolution density 
field at a redshift z = 0, for the model L256 using the parameter set A 
(see Sect. 2.5 for explanation); 

Column 4: The number of DF clusters in the high-resolution density 
field at a redshift z = 0, for the model L256 using the parameter set B; 
the number of AHF halos at the redshift z = for the models M256 and 
LIOO. 

All models of the same series have the same realisation of 
random Fourier amplitudes, so the role of difl'erent waves in the 
models can be easily compared. To compare difl'erent models of 
a series, every particle has an identification number, the same 
for all models of the same series. The main model parameters 
are given in Table [TJ In addition to models with cuts on large 
scales, we calculated a model, where the density perturbations 
on intermediate scales were truncated to zero at the wavelengths 
8-64 hr^ Mpc, the model M256.864. This model shows the 
eff'ect of the absence of perturbations of medium scales. 

We use the following notation for our models: the first char- 
acters M and L designate models with resolutions of A^grid = 256 
and A^grid = 512, respectively; the following number gives the 
size of the simulation box, L, in hT^ Mpc; the subsequent num- 
ber indicates the maximum wavelength used in the simulation, 
also in hT^ Mpc; if the full power spectrum is used these two 
numbers coincide. The locations of the cells inside the cubical 
density grid are marked by cell indices (/, 7, K), where /, j, and k 
are integers that run from 1 to A^grid- 
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M256.256 z = 30 M256.256 z = M256.864 




Fig. 2. The left and middle panels show the power spectra for the models of the series M256 at the epochs z = 30 and z = 0.00, 
respectively. The right panel shows the power spectra for the model M256.864 at the epochs z = 30 and z = 0.00. In the middle 
panel, we also show the power spectrum for the full model M256.256 at the present epoch z = where the evolution has been linear 
on all scales. 



For the models of the M256 and the M768 series, we 
used in simulations the AMIGA code, which is the follow-up 
of the multi-level adaptive particle mesh (MLAPM) code by 
iKnebe et aP ("2001). This code uses an adaptive mesh technique 
in the regions where the density exceeds a fixed threshold. In 
this code, gravity is automatically softened adaptively, so that 
the softening length is close to its optimum value in both high 
and low-density regions. We chose a maximum level of eight 
refinements. For the models of the LI 00 and L256 series, we 
used the GADGET-2 code with gravitati onal softening length 
10 kpc and 20 kps, respectively (ISpringel et al.L 
ISpringelL 1200^ . The simulations M256 and M768 were per- 
formed at the Tartu Observatory, the simulation LI 00 at the 
Astrophysical Institute Potsdam, and the simulation L256 at the 
High Performance Computing Centre of the University of Tartu. 

The initial density fluctuat ion spectra w ere generated using 
the COSMICS code by Bertschinger| (ll995lPl . We assumed the 
cosmological parameters Qm = 0.28, Qa = 0.72, erg = 0.84, 
and the dimensionless Hubble constant h = 0.73; to generate the 
initial data, we used the baryonic matter density Db = 0.044 
(iTegmark et al.l (l2004l) ). We generated initial conditions for a 
truncated power spectrum, using the full spectrum to calculate 
the model parameters. The amplitude of a spectrum was set to 
zero for k < kt during the calculation of the initial density field, 
keeping all simulation parameters fixed across the full set of 
realisations. The calculations were started for an early epoch, 
z = 100. The particle positions and velocities were extracted for 
12 epochs between the redshifts z = 30 ... (the models of the 
M256 and the M768 series), and for 7 epochs between the red- 
shifts z = 100 ... (the models of the LlOO and the L256 series). 

2.3. Calculation of the density field 

For every particle, we calculated the local density in units of the 
mean density, using the positions of 27 nearby particles. We also 
calculated the mean global density at the location of the particle. 
For this purp ose, we first determined the density field, using a 
B3 spline fsee lMartinez & Saail.[20Q2b : 

B3(x) = ^ - 2P - 4|x - Ip + 6\x\' - 4\x + l\' + |x + 2\'] , (4) 

where function diff'ers from zero only in the interval x e [-2, 2]. 
The one-dimensional B3 box spline kernel of width h = N is 

K^^\x;N) = B3(x/N)/N. (5) 

^ |http : //arcturus .mit . edu/cosmicsl 



This kernel preserves the interpolation property (mass conserva- 
tion) for all kernel widths that are integer multiples of the grid 
step, h = N. The 3D K^^^ box spline kernel is given by the direct 
product of three one-dimensional kernels 

Kb(x; N) ^ Kf(x; N) = 4'^)(x; N)K^^\y; N)K^^\z; N), (6) 

where x = {x,y,z}. 

To calculate the high-resolution density field, we use the 
kernel of the scale equal to the cell size of the simulation, 
h = L/Ngrid, where L is the size of the simulation box, and 
A^grid is the number of grid elements in one coordinate. We use 
consecutive smoothing stages; the smoothing stage of the or- 
der / has the smoothing scale n = L/Ngnd x 2^ The eff'ective 
scale of smoothing is equal to 2 x r/. We applied this smooth- 
ing up to index (order) 6. The smoothing with the B3 spline 
kernel is rather close to the smoothing with an Epanechnikov 
kernel of the same scale, ii sed in earlier studies of superclusters 
(lEinasto et al.Ll2006ll2QQ7h . For the model LlOO, the smoothing 
index 5 corresponds to the kernel of the radius 6.25 Mpc, for 
the model L256 the smoothing index 4 corresponds to the kernel 
of the radius 8 Mpc, and for the model M768 the smooth- 
ing index 2 corresponds to the kernel of the radius 12 Mpc. 
These kernel radii were used to calculate the global density field 
used in the selection of the void and supercluster core regions. 

The power spectra of the models of the series M256 are 
shown in Fig.[2]for both an early epoch, z = 30, and the present 
epoch z = 0.0. The power spectra for the model M256.864 
are shown in the right panel of the same figure for the same 
epochs. For comparison, we also show in the middle panel the 
power spectrum of the model M256.256 for the present epoch, 
if evolved linearly at all scales. Fig. [2] shows that at an early 
epoch and at small wavelengths all spectra practically coincide. 
At the present epoch, the amplitudes of spectra for diff'erent 
models diff'er at small wavelengths: the models with the trun- 
cated spectra evolve at small wavelengths more slowly than the 
full model. At small wavelengths, the slowest evolution is seen 
in the model M256.008, where the growth of the amplitudes is 
almost linear on all scales. This eff'ect is expected because den- 
sity waves of larger size that could amplify the growth of small- 
scale waves are absent (see the analysis below). We also see that 
in all models with a cutoff', the amplitude of the spectrum in 
the cutoff' range slowly increases, i.e. some power is translated 
from smaller to larger waves. This eff'ect is much stronger in the 
model M256.864, where in the cutoff' range of the spectrum a 
fairly high amplitude of the spectrum develops for the present 
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Fig. 3. Density fields for the models of the M256 series. The upper panels show the high-resolution fields for the models M256.256, 
M256.064, and M256.032, the lower panels - for the models M256.864, M256.016, and M256.008 (from left to right). The densities 
are shown for a layer of 6 /z"^ Mpc thickness at the k = 15 coordinate. All fields correspond to the present epoch z = 0. The densities 
are expressed on a logarithmic scale to enhance the low-density regions. Only the overdensity regions are shown. The upper limits 
to the density are 40 in the upper panels, and 10, 10, and 5 in the lower panels. 



epoch, only a factor of ten lower than the amplitude at the cutoff' 
wavelength A = 8 /z"^ Mpc. 

The high-resolution density fields at the present epoch for 
models of the M256 series are shown in Fig.[3l 

2.4. Comparison of halo finders 

In simulations, instead of groups and clusters of galaxies, dark 
matter halos are usually used to define simulated compact sys- 
tems. In the pilot phase of this study, we used the PM code in 
a box of size L = 128 /z"^ Mpc with 128^ particles and simu- 
lation cells. For these models, we applied the conventional halo 
finder using the FoF technique with the linking length parame- 
ter 0.2 in units of the mean separation of particles. The analysis 
showed that the FoF halos are of little use for our study, since in 
the central regions of superclusters, the richest FoF halos contain 
the whole core of the supercluster, and have sizes of the order of 
10 r 1 Mpc. 

For this reason, we tri ed the adaptive Amiga Halo Finder 
(AHF) code developed bv iKnollmann & Knebel ('2009'). where 
the number of particles in a halo is A^p > 20. Halos and their 
parameters (masses, virial radii, positions, velocities etc.) were 
found for most models and simulation epochs. Here we have an- 
other difficulty: the number of halos for early epochs (high red- 
shifts) decreases dramatically with increasing z: there are about 
1000 times fewer halos at the redshift z = 10 in the model LI 00 
than at the present epoch. For models of the M256 series, the 
AHF ffnds no halos at the redshifts z = 5, and 10, about 10 times 



fewer halos than the density ffeld cluster method (see below) for 
the models M256.256, M256.064, M256.032, and M256.016, 
and almost no halos for the model M256.008 (see Table [B. The 
density enhancements in these models are much smaller and do 
not conform to the deffnition of the AHF halos, (which have at 
least 20 particles), in spite these models containing many small 
compact density enhancements at all redshifts. In Fig. [S] we 
present the density ffelds of the models of the M256 series at the 
redshift z = 0, and Figs. 1-3 of Einasto et al. (2010) for higher 
redshifts. 

Taking these difficulties into account, we used the high- 
resolution density field to define compact systems, i.e., den- 
sity field (DF) clusters. DF clusters correspond to groups and 
clusters of galaxies in the real Universe, and to halos in sim- 
ulations. In defining density field clusters, we applied the con- 
vention al trad i tion, used by Abe l l in h is rich cluster cata- 
logues (lAbellL 119581: lAbelletaLl Il989h . The Abell clusters 
were selected using galaxies in a sphere of a radius of about 
1.5 h'^Mpc. Our density field m ethod to find halos is rather 
simi lar to the DENMAX method by lBertschinger & Gelbl(ll991h 
and Gelb & Be rtschinger ('1994'). The basic diff'erence lies in the 
detail: we do not use particles to define the center and the mass 
of the halo, but only data obtained directly from the density field. 
The procedure for finding the DF clusters is described in the next 
subsection. 

To compare our DF cluster mass distributions with the mass 
distributions of the AHF halos, we applied both methods for 
the highest resolution model LI 00. We found the DF clusters 
using the parameters Dq = 2 (in units of the mean density). 
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Dp = 20 (the minimal number of particles), and N = 8; the 
latter corresponds to the virial (outer) radius of the halo Ryi^ = 
1.66 Mpc, which is similar to the virial radii of the largest 
halos found with the AHF method. The cumulative mass func- 
tions for the DF clusters and the AHF halos for all models of the 
LI 00 series for the redshift z = are shown in Fig.lH 

Fig.|4]shows that the maximum masses of halos for both halo 
finders are practically identical for all models with diff'erent cut- 
off's. The shapes of the cumulative mass functions are slightly 
different. Table[T]and Fig.|4]show that the total numbers of halos 
found with the AHF method are about three times larger than 
the numbers of the DF clusters for models with the same cut- 
off. This difference is expected, since AHF treats all subhalos as 
independent halos. 

In the mean mass range, the number of the DF clusters is a bit 
higher than the number of the AHF halos. This difference is also 
expected. Our DF cluster finder always has an identical volume 
around a high-density peak, whereas the virial radii of the AHF 
halos decrease for the model LlOO.lOO from 7?vir = 1.5 Mpc 
to /?vir = 0.05 Mpc for the considered mass range. For this 
reason, close intermediate-mass halos are counted as separate 
small halos in the AHF method, but as a single larger halo in the 
DF cluster method. 

These differences in the mass distributions are not important 
for the present study, since we compare mass distributions for 
different cutoffs of the power spectra, obtained with the same 
halo identification method. We use halos basically to identify 
voids and to measure their radii; the internal structures of the 
halos are not important for the present analysis. 

2.5. Density field clusters 

As a first step, we find local density maxima of the high- 
resolution density field. In the search for maxima, we use three 
parameters: the minimum density threshold Dq, above which we 
search for density maxima, and the minimum mass of the DF 
clusters. Dp. During the search, densities are expressed in units 
of the mean density of the simulation, and masses - in units of 
the particle mass of the simulation (the number of particles in a 
halo). The mass of the DF clusters, Md, is calculated by adding 
the local densities in the cells within ±N cells from the central 
one, i.e. in a total in (2N -h 1)^ cells. The number N defines the 
volume where we count the mass of the cluster, and it is the third 
parameter of the search. We express the DF cluster masses in so- 
lar units, using the known particle mass. 

For the M256 model, we used the parameters Dq = 2, 
Dp = 5, and A/^ = 3; for the M768 model, we used the param- 
eter set Do = 2, Dp = 5, and N = I, which correspond to the 
cluster search radii r = 3.5 Mpc and r = 4.5 Mpc, re- 
spectively. To find the dependence of the DF cluster mass func- 
tions and the cluster-defined void radii on the parameters used 
in the definition of a DF cluster, we used two sets of parame- 
ters for the models of the L256 series - a high- sensitivity set A 
with Do = 2, Dp = 5, and N = 3, and a low- sensitivity set B 
with Do = 10, Dp = 50, and N = 5. The cluster search radii are 
r = 1.75 Mpc, and r = 2.75 Mpc, for the parameter sets 
A and B, respectively. The parameter set A finds the DF clus- 
ters in faint filaments crossing supercluster-defined large voids, 
the parameter set B avoids most of the DF clusters in these fila- 
ments. 

The total numbers of the DF clusters found for the present 
epoch z = and the selection parameters used are given in 



Table [T] For the models M256 and LI 00, we also give the num- 
bers of the AHF halos. 



3. Analysis of models 

3.1. The density field of truncated models 

To have an idea of the appearance of the density fields, we show 
in Fig. [3] the high-resolution density fields for all the models of 
the series M256: M256.256, M256.064, M256.032, M256.016, 
M256.008, and M256.864. To emphasise the filaments joining 
the DF clusters, we express densities in logarithmic scale, and 
use sheets of a thickness 6 Mpc. The fields are shown for the 
present epoch z = 0.0. The models M768 and L256 look similar, 
when cuts at the same scales are used, but fewer details are seen 
in the model M768 with respect to the model L256. 

Fig. [3] shows that the models M256.256 and M256.064 dis- 
play rather similar patterns of the cosmic web. An even closer 
similarity is observed when we compare the density fields of the 
models L256.256 and L256.128. This similarity of the patterns 
of the cosmic web in these models shows that the pattern of the 
web is defined by density perturbations of scales smaller than 
^ 100 h-^ Mpc. 

In models with shorter scale power spectrum cuts large- 
scale features disappear. In the models with a cutoff scale < 
100 Mpc, the shape of the pattern of the cosmic web is de- 
termined by the maximum scale of density perturbations. This 
is one of our main qualitative conclusions: in the absence of 
large-scale perturbations, galaxy systems larger than the cut- 
off scale do not form. In other words, the scale of the cosmic 
web is determined by the density waves of the largest scale 
present. However, this conclusion is correct only up to the scales 
~ 100 Mpc. Fig.[3](and similar plots for other models) shows 
that the addition of perturbations larger than this scale does not 
create larger systems, but only amplifies them. We study this 
phenomenon below. 

3.2. The role of the density waves of medium scales 

To see the effect of the absence of density perturbations of 
medium scales, we performed one simulation, where the am- 
plitude of the initial density fluctuations between wavelengths 
8 /z"^ Mpc and 64 Mpc was truncated to zero. This model 
(M256.864) is otherwise identical to the model M256.256 and 
was generated using identical Fourier amplitudes. The power 
spectra of this model at z = 30 and z = are shown in Fig. [2l 
and the density field at the present epoch in Fig.O 

Fig. [3] demonstrates very clearly the role of density perturba- 
tions of various scales in the formation of the cosmic web. In the 
absence of large-scale perturbations, systems larger than the cut- 
off scale do not form. The result of removing only the medium- 
scale perturbations is of particular interest. As seen from Fig.O 
during the evolution the power spectrum on intermediate scales, 
which is absent at the initial epoch, is in this case almost re- 
stored. However, this increase in the power spectrum amplitudes 
does not lead to the formation of medium-scale galaxy systems 
such as filaments. In a model without the initial medium-scale 
perturbations, filaments are absent, both within superclusters and 
between them. The distribution of small-scale systems is more or 
less random, and there are no compact systems of galaxies such 
as clusters - the compact systems are rather small. 
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LlOO LlOO 




Fig. 4. The left panel shows the cumulative mass functions of the density field clusters for the model LlOO for various cutoff' scales. 
The right panel shows the cumulative mass functions of the AHF halos for the same model. 




Fig. 5. The matrix of distances to the nearest DF clusters for the models M256.008, M256.016, and M256.256 (from left to right). 
The "density" is defined as the distance to the nearest DF cluster. The matrix is given for the k = 220 Mpc coordinate. Note 
that changes of the mean radii in the distance matrix are similar to changes of the overdensity pattern of the density field, shown in 
Fig.S 



3.3. The distribution of void sizes 

Cosmic voids are regions of space devoid of certain kinds of 
objects - galaxies, clusters of galaxies etc. Diff'erent types of 
objects define voids of diff'erent size. The reason for the de- 
pendence of void sizes on the mass (or the luminosity) of ob- 
jects used in their definition is simple. Large voids are deter- 
mined by rich clusters and crossed by filaments of faint galax- 
ies. Moreover, almost all systems of galaxies contain outlying 
faint members (see Fig. [T] for the luminosity density field of a 
spherical shell of the Sloan Digital Sky Sur vey). Dwarf galaxies 
define much smaller voids t han g iant on es (lEinasto et al.L 1986 . 
19891, '1991'; 'Lindner et al.', 1995, ^199^ iGottlober et al.L [2003 : 
von Benda-Beckmann & Miiller, 200 8). For the hiera r chy o f 
voids, ^y e also r e fer to Ivan de Wevgaert & v an Kampen ( 1993), 
iPeeblesI (|200ll), IGottlober et al.l (<2003l). 'Arag on-Calvo et al 
(2007), van de Wevgaert & Platen (2009), Ara gon-Calvo et al.l 
(i2010al) . and lAragon-Calvo et al.i 62010b) . 

To find the distribution of void radii, we used a simple void 
finder suggested by Einasto et al. (1989). For comparison of dif- 
ferent void finders, we refer to IColberg et al.l (120081) . For each 
vertex of the simulation grid, we first calculated its distance 
to the nearest DF cluster. For the positions of the DF clusters, 
we used the /, 7, ^-indices of the maximum local density cells. 
The distance matrix is similar to the density field matrix; for the 
models M256.008, M256.016, and M256.256, these matrices are 
shown in Fig.[5j 



Fig. [5] shows the maxima of the distance matrix field. These 
maxima correspond to the centres of voids, and their values are 
the void radii. The distribution of clusters is noisy, thus there 
are many nearby local maxima in the distance matrix field. We 
define the position of the void centre as the location of the cell, 
which has the highest distance value within a box of the size of 
±3 grid elements. 

Fig. [5] shows that in some places there are long ridges in the 
distance matrix; these ridges delineate elongated voids. We con- 
sider them to be individual voids, if their centres are separated 
by more than 3 Mpc, i.e. they count as separate entries in the 
void search using the criterion shown above. 

In models with a higher cutoff', the maximum void radii in- 
crease with the cutoff' scale, but only moderately. There is only 
a small diff'erence in the distribution of void radii between the 
models M256.256 and M256.064. In other words, very large 
density waves do not change the structure of the supercluster- 
void network, their role consisting essentially of the amplifica- 
tion of existing structures. 

The mean void radii were found for a broad range of the DF 
cluster mass thresholds, Mth. The highest threshold used was se- 
lected so that the volume density of the DF clusters in the sample 
is approximately equal to the volume de nsity of the Abel l clus- 
ters, at about 25 x 10"^ (/z"^ Mpc)^ (Ein asto etall l2006h . This 
gives ^ 420 DF clusters in the volume of the L256 model. The 
actual number of DF clusters used in the void definition depends 
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Mhalo [10^^12 Msun/h] 



Mth [lO'^lOMsun/h] 



Fig. 6. The left panels show the cumulative mass functions of the density field clusters for the models with various cutoff' scales. The 
right panels show the mean radii of voids, defined by the DF clusters for diff'erent threshold masses, Mth, and for various cut scales. 
The upper panels are for the model L256 with the DF cluster search parameter set A, the middle panels are for the same model with 
the DF cluster search parameter set B, the lower panels are for the model M768. 



on the value of the DF cluster mass threshold, Mth, used in the 
void search. 

In Fig. [6l we show the mean void radii as functions of the 
threshold mass Mth of the DF clusters used in void search. This 
figure shows that the void radii depend on the threshold mass, as 
expected from similar studies of the void radii in galaxy samples 
of diff'erent threshold luminosity, demonstrated in many studies 
cited above. We study this eff'ect in more detail below. 

3.4. The distribution of the DF cluster masses 

One of the characteristics of the cosmic web is the distribution of 
cluster masses. The cumulative mass functions of the DF clusters 
for the models of the series L256 and M768 are shown in Fig.[6j 
We see that DF cluster masses strongly depend on the scale of 
the power spectrum cutoff'. All models with a cutoff' on the scale 



8 Mpc have DF cluster mass distributions with rather sharp 
decreases on the high mass side. The maximum masses of the 
DF clusters in these models are ^ 1.5 x 10^^ M© for the L256 
model and ^ 4 x 10^^ M© for the M768 model. The maximum 
masses of the DF clusters in the models of the LI 00 series are 
lower than the maximum masses of the DF clusters in the models 
of the L256 series. These diff'erences are due to the larger sizes of 
the boxes used in the DF cluster search in the models of the L256 
and the M768 series. As shown in Fig.|4l there are practically no 
diff'erences between the maximum masses of the DF clusters and 
those the AHF halos of the LI 00 model, since in the DF cluster 
search we used a search radius approximately equal to the virial 
radius of the most massive AHF halo. 

With increasing spectrum cutoff' scale, the mass distribu- 
tions rapidly shift to higher masses. This rapid increase in the 
maximum DF cluster mass continues up to the cutoff' scale of 
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64 Mpc. A higher cutoff scale for the power spectrum only 
moderately increases the maximum masses of the DF clusters. 
This increase almost stops at the cutoff scale 128 Mpc. 
The most massive DF clusters in these models have the masses 
^ 7 X 10^4 and ^ 2 X 10^^ M©, for the models L256 and 
M768, respectively. 

The number of low mass DF clusters in the models depends 
on the selection parameters used in the cluster search. The total 
number of DF clusters in the model L256 found for the high- 
sensitive and the low-sensitive parameter sets, A and B, is rather 
different. The high- sensitive set A with a lower search peak den- 
sity level Dp = 5 has about six times more DF clusters than the 
set B with Dp = 50. All additional DF clusters found for the 
parameter set A are in the low mass range, as seen in Table [T] 
and Fig.O In the higher mass region, the cumulative mass func- 
tions for both sets of search parameters almost coincide. As we 
see below, in the void search both parameter sets yield voids of 
approximately equal mean radii for identical threshold masses. 
This shows that our results are rather robust to the parameter 
choice in the DF cluster definition. 

A more detailed study shows, however, that density waves of 
still larger scales influence the structure of the cosmic web. We 
discuss below the evolution of the density field in the void and 
supercluster core regions. 

4. Discussion 

4.1. Voids defined by clusters of different mass 

As seen in Table [T] and Fig. [6l the numbers of the DF clusters 
for the search parameter set A for the models L256 are about six 
times larger than for the parameter set B. The basic reason for 
this difference is the use of a much lower threshold density Dq = 
2 and the peak density Dp = Sin the cluster search; for the set B, 
the corresponding densities are Dq = 10, Dp = 50. For instance, 
the number of the DF clusters in the L256.256 model for the 
lowest mass threshold 10^^ M© is 149265 for the parameter set 
A, and 59100 for set B. In spite of the large difference in the 
number of the DF clusters for different search parameters, the 
mean void radii for both parameter sets are approximately equal, 
as seen in the right panels of Fig.[6l A detailed inspection of the 
density field maps shows that this similarity in void sizes has 
a simple reason: the DF clusters of the set B are high-density 
peaks inside the same filaments that form the basic web of more 
densely populated filaments of the DF clusters of the set A. The 
density field maps also show that there are no extra faint clusters 
inside the relatively small voids defined by the clusters of the 
set A. For void definition, it is irrelevant how densely populated 
the filaments surrounding the voids are, if there are no galaxies 
outside the filaments. 

This conclusion agrees with the results of earlier stu dies by 
iLindner et al.l(ll996h JPeebles (2001), and lGottlober et al.l(l2003l) 
among others, that voids defined by relatively faint galaxy fila- 
ments are completely devoid of galaxies. 

Fig. [6l shows that for a low DF cluster mass threshold, the 
void radii are almost independent of the mass threshold Mth. 
This means that in this DF cluster mass interval the clusters are 
located in identical filaments. 

If the DF cluster mass threshold increases further, the void 
radii start to increase. This means that some filaments are fainter 
than the respective mass threshold limit, and do not contribute to 
the void definition. This slow increase in the mean void radii Ry 
as a function of the DF cluster mass threshold for strongly 
cut models is valid until a certain value, of about 3x10^^ M© 



and 10 X 10^^ M© for the L256.008 and the L256.016 models, 
respectively. Thereafter, with an increasing mass threshold, the 
void radii increase very rapidly until clusters disappear. This ef- 
fect is due to the very sharp decrease in the number of DF clus- 
ters of a high mass (see the left panels of Fig. [6]). These rare 
clusters define very large voids. The sizes of these voids are not 
characteristic of the overall cosmic web pattern of the particular 
model. 

Fig. [6] shows that in the models with a higher cutoff scale 
Acut > 32 Mpc, the void mean radii grow with the growth 
of the mass threshold smoothly, until the highest masses, which 
correspond to very rich clusters. Moreover, the Ry versus 
curves for the models M256.032, M256.064, M256.128, and 
M256.256 are almost identical. This means that the scale of the 
cosmic web is determined essentially by density perturbations of 
a scale up to 32 /z"^ Mpc. Some differences in the Ry versus 
curves remain between the models M256.032 and M256.064. 
The higher cutoff models are practically identical in this void 
size test. 

Thus, the void analysis confirms our results from the mass 
distribution of the DF clusters, that density perturbations of large 
scales have little effect on the pattern of the cosmic web as char- 
acterised by void sizes. 

The largest voids in the L256.256 model are defined by the 
rich clusters of mass threshold 9x10^^ M©, and have the mean 
radii ^ 30 Mpc. These voids are characteristic of filaments, 
which divide supercluster-defined voids (supervoids) into ap- 
proximately equal sub- voids. 

The model M768 has a lower resolution; here faint filaments 
are absent, and the scale of voids is defined by the DF clusters of 
higher mass. As in the model L256, models with a smaller cut- 
off scale have smaller voids. The mean void radii for the models 
M768.128 and M768.768 have practically identical Ry versus 
Mth curves. The largest voids in this model are defined by the 
DF clusters of the mass threshold 6 x 10^"^ M©. These very rich 
clusters populate the cores of rich superclusters, thus the respec- 
tive void radii ar e practically equal to the radii of supervo ids, 
Ry ^ 10 h-^ Mpc (I Lindner et allll995l:lEinasto et al.Lll997ah . 

4.2. The evolution of the density field in the void and 
supercluster core regions 

As an additional test, we studied the evolution of the density 
field in extreme voids and supercluster cores. Here we used the 
property of our models - all variants of a series with a different 
spectrum cutoff scale were generated with an identical random 
number set. Thus we can use the particle ID numbers to study 
the behaviour of the web on different spectrum cutoff scales. 

To define the extreme void and supercluster core regions we 
used the global density field smoothed with a wider kernel, ap- 
plying the same procedure as explained above for finding the 
DF clusters in the high-resolution density field. To calculate the 
global density field, we used for the models of the series L256 
the kernel size 8 /z"^ Mpc, for the models of the series M768 the 
kernel size 12 Mpc when applying the B3 spline. For every 
particle, we store the local density value at the particle location 
(found on the basis of 27 nearest neighbours) in addition to its 
coordinates, and the global density found from the density field 
with a large smoothing length as described above. These data are 
available for all models and evolution steps. We also found the 
distributions of the local and global density of particles. 

Using this information, we extracted particle samples for 
models with the full power spectra at the present epoch, which 
had ^ 5 % of particles with the lowest and the highest global 
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Fig. 7. The cumulative distributions of the local densities for particles in the void and core regions for the present epoch z = 
are given in the left and the right panels, respectively. The upper row is for the L256 models, the lower row for the M768 models. 
Models with various cutoff scale of the power spectra are shown. 



density values. We refer to these particle samples as the void and 
core samples, respectively. The number of particles in the void 
and core samples of the M768.768 and the L256.256 models are 
given in Table [2j We calculated the cumulative distributions of 
the local densities of the void and core particles at the present 
epoch z = for all the models of the L256 and M768 series. The 
results are shown in Fig.[71 

We discuss first the cumulative density distribution of the 
void and core particle samples of the M768 model. This model 
has a lower resolution and lacks the low-mass DF clusters in the 
deep void regions, so the interpretation of the results is simpler. 
Fig. [7] shows that the cumulative particle density distributions 
of the void and core regions of the model M768.012 are almost 
identical. Both distributions are rather symmetrical around the 
mean density level Aoc = 1 • This result shows that in this model, 
there is no diff'erence between the void and core regions. The 
symmetry around the mean density level indicates that the evo- 
lution is still close to the linear regime of the growth of density 
perturbations, i.e the growth of the negative and positive sections 
of the density contrast S = D - I is similar and proportional to 
the growth factor of the evolution. 

In the model M768.032, the amplification of the density 
growth by the density waves of medium scale is already visi- 
ble. The diff'erence between the density distributions of the void 
and core samples is still small - there are no big voids and no 
rich supercluster cores. 

The distributions of the particle densities in the void and core 
regions are completely diff'erent in the models M768.128 and 
M768.768. In the void regions, all local densities are lower than 
the mean density, i.e. there are no systems of galaxies. We recall 



that for the formation of a galaxy or a galaxy group the local 
density must exceed a certain threshold, about 1.6 in the mean 
density units. In contrast, in the core regions there are almost no 
particles of local densities less than the mean density. 

Table 2. Results of the K-S test. 



Model n d p 



M768 void 820002 0.07563 0.000000 

M768 core 823717 0.12304 0.000000 

L256void 6677929 0.01747 0.000000 

L256core 6768219 0.04725 0.000000 



The meaning of n, d and p is explained in text. 

We also note that there is a small but definitely significant 
diff'erence between the distribution of densities in the void re- 
gions of the models M768.128 and M768.768, and in the core 
regions of the same models. The diff'erence is smaller in the void 
regions, but in both regions it is clearly present. The results of the 
Kolmogorov- Smirnoff' test of the comparison of the cumulative 
distributions of particle densities in the void (core) regions of the 
M768.128 versus M768.768 and the L256.256 versus L256.128 
models are shown in Table |2l where n is the number of parti- 
cles in respective samples (which are identical in the full model 
and the 128 Mpc cut models of the same series), d is the 
maximum diff'erence of cumulative distributions, and p is the 
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probability that the distributions compared are taken from the 
same parent sample. As we see, this test shows that the distribu- 
tions are different at high significance levels. This shows that the 
density perturbations of a larger scale than 128 Mpc make 
voids emptier and systems of galaxies richer, i.e. they amplify 
the emptiness and the richness of the cosmic web. 

The model L256 has much higher spatial and mass resolu- 
tions. For this reason, there are more diff'erences between the 
particle local density distributions of the void and core regions. 
In this model in the void regions we see the presence of particle 
overdensities, with local densities D > 1. These particles form 
faint filaments crossing large voids. It is important that the larger 
the power spectrum cutoff' scale (i.e. the larger perturbations are 
included in the simulation), the lower the fraction of the number 
of particles with overdensities. In other words, large-scale per- 
turbations suppress the formation of filaments in void regions. 

In the core regions the growth of overdensities is more rapid, 
the larger is the scale of the perturbations included in simula- 
tions. In the model L256.256, there are no particles with local 
densities below the mean density, i.e. all particles have been 
moved to rich systems. Owing to the higher resolution, the high- 
est densities in the core regions are much higher than in the 
model M768. 

4.3. Why the perturbations of the largest scales do not 
increase the scale of the cosmic web? 

Our study suggests that density perturbations up to the scale 
^ 100 /z"^ Mpc determine the scale of the cosmic web in terms 
of void sizes. In contrast, waves of larger wavelengths do not in- 
ffuence the scale of but only amplify the web, leaving its scale 
unaff'ected. Thus, it is unclear, why the growth of the scale of the 
skeleton of the cosmic web with increasing wavelengths of den- 
sity perturbations, stops at the perturbation scale ^ 100 /z"^ Mpc. 
Has this some deeper physical meaning? That the change in the 
behaviour of the waves of diff'erent scale occurs in numerical 
simulations where only standard physics of the early Universe is 
included, implies that there is a simple answer. 

For h ^ 0.7, this largest pattern scale is close to 7?eq = 

ci(to)r]Qq = 2( V2-l)(c/ifo)x V^rad/^m ~ 120 Mpc, which is the 
only scale characterising a primordial scale-free (^s ~ 1) spec- 
trum of density perturbations multiplied by the scale-dependent 
transfer function arising due to the transition from the radiation- 
dominated stage to the matter-dominated one, which occurred 
at z = Zeq ~ 3200 according to the most recent observa- 
tional data (see [Jarosik et al.l (|2010)). Here the conformal time 
dz/H(z), H = a(t)/a(t) is the Hubble pa- 
rameter, to is the present time (H(to) = Hq), and Drad includes 
the contribution from three species of neutrino which all may be 
considered massless at the moment of matter-radiation energy- 
density equality z = Zeq. We note that Rq^ does not depend on h 
if expressed in terms of Qradh^ ^ and ^mh^, nor depend on 
the dark energy equation of state wde (^de = - 1 for a cosmo- 
logical constant), only the present value of Qm being important. 
Moreover, the expression for /?eq remains the same for an open 
(negatively spatially curved) Universe in the absence of a cos- 
mological constant. 

Additional insight to this problem is given by simulations 
designed to understand the possible future of the develop- 
ment of the cosmic web . These s imulations were performed 
by a number of authors (Loe bl 12002 : Nasamine & Loeb, 2003; 
Bushaetal, 2005; Diinner et al.L I2006L "120071: iHoffman et all 
i2007.:.Krauss & Scherrer^ .20071) . In the early Universe, it is well 



known that the matter density Dm dominated over the energy 
density. As time proceeded and the Universe expanded, the Qa 
term increased and caused the Universe to expand in an acceler- 
ating fashion from a redshift of z ^ 0.8. 

Our results for the absence of the growth of the scale of the 
cosmic web with the increase in the scale of density perturba- 
tions over ^ 100 Mpc can probably be explained as the re- 
sult of the freezing of the web that started at recent redshifts, 
I -\- z < ((1 - ^lm)/^m)^^^ for an exact cosmological constant. 
Simulations of the future development of the web mentioned 
above support this interpretation. As stated by iHoff'man et al.l 
( 2007) "in comoving coordinates the future large-scale structure 
will look like a sharpened image of the present structure: the 
skeleton of the cosmic web will remain the same, but clusters 
will be more isolated and the filaments w ill become thin ner". 
Simulations by Naga mine & Loebl (l2003l) and floff'man e t al.l 
(2007) among others have shown that the mass evolution of 
bound objects such as clusters will stop. In the very distant fu- 
ture, all galaxies that are not bound to the Local Group (which 
consists basically of the merged Milky Way and M31) will fade 
away while approaching the event horizon, as seen from an ob- 
server inside the Local Group. It remains to be investigated how 
a possible deviation of wde from -1, if it exists, can aff'ect the 
freezing and the properties of the cosmic web. 

5. Conclusions 

The basic conclusions of our study are as follows: 

- The properties of the cosmic web depend strongly on density 
perturbations of various scales. 

- Small-scale perturbations up to the scale ^ ^ Mpc are 
responsible for the formation of galaxy and cluster type sys- 
tems. 

- Medium-scale perturbations of the scale 8... 64 h'^Mpc 
form the filamentary web inside and between the superclus- 
ters. They also contribute to the amplification of the systems 
formed by perturbations of smaller scales. 

- The cosmic web with filamentary superclusters and voids is 
formed by the combined action of all perturbations up to the 
scale ^ 100 Mpc. The largest perturbations in this range 
determine the scale of the supercluster-void network. 

- Perturbations of the largest scales > 100 Mpc modulate 
the richness of galaxy systems from clusters to superclusters, 
and make voids emptier. 
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